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We study low energy spin and charge transport through a ID Mott-Hubbard insulator of finite 
length L attached to Fermi liquid reservoirs characterized by different chemical potentials for elec- 
trons of opposite spin polarizations as it happens in quantum spin Elall insulators. We calculate the 
average currents (charge and spin) and their correlators and demonstrate how a transition induced 
by the reservoirs to the low energy Fermi liquid regime results in breakdown of the spin-charge 
separation, which is visible in the presence of the spin dependent voltages and a weak one electron 
scattering in the system. These calculations are carried out under assumption that the Hubbard 
gap 2M is large enough: M > Tl = Vc/L {vc- charge velocity in the wire) and the scattering rate 
Vs <^Tl- Relation of these results to Kondo dot transport in the Toulouse limit is also clarified. 
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Spin-charge separation in a ID Mott-Hubbard insu- 
lator (MHI) known from the exact solution of the ID 
Hubbard model [l| at half-filling and confirmed in exper- 
iments with quasi ID materials 0, 0] can lead to an un- 
usual effect if the insulator is used for transport between 
Fermi liquid (FL) reservoirs filled with electrons of oppo- 
site spin polarization up to different chemical potentials. 
Different spin-dependent chemical potentials can arise, in 
particular, in 2D quantum spin Hall insulators (QSHIs) 

, where transport is carried by pairs of edge states of 
up and down spin polarized electrons moving in opposite 
directions. Since in MHI the charge transfer realized by 
holons carrying one-electron charge e (e = ft = 1 below) 
is suppressed by a Hubbard gap 2M, but spin carrying 
excitations, spinous, move freely, the perfect MHI serves 
as an ideal spin current transmitter. Any small difference 
in the spin dependent voltages V^^a = ± reverses direc- 
tion of one of the spin polarized currents in spite of large 
average difference between the source and drain chemi- 
cal potentials. In the case of the QSHI reservoir the spin 
polarized tunneling currents can be found by measuring 
edge currents to the right and left from the juncture with 
the MHI wire along the reservoir boundary. If there is 
only a single pair of edge states in the reservoir and the 
spin dependent voltages are applied anti-symmetrically, 
the ideal spin current transmittance through the MHI 
makes both these edge currents and, hence, conductance 
along the reservoir boundary zero, since the difference in 
the currents of the outgoing spin polarized edges is op- 
posite to the difference of the incoming currents. This 
simple picture of the transport can be spoilt, however, 
by one electron processes which affect spin-charge sepa- 
ration in the MHI. 

In this paper we consider effect of the FL reservoirs on 
low energy spin-charge separation in the MHI of finite 
length L: M > Tl = Vc/L {vc- charge velocity in the 
wire) in the presence of one-electron impurity backscat- 
tering of low rate Fs ^ M. For low energy (< Tl) the 



problem has been mapped [sj through a Duality Trans- 
form onto the model of a point scatterer with pseudospin 
imbedded in TLL and solved by fermionization. This 
model is also related to the Toulouse limit in Kondo dot 
transport @ as discussed below. Its solution shows a 
Kondo- like resonance of conductance below the crossover 
energy F = Fj -I- Fc, which remains exponentially small 
(Fc « y/TlM eiq){-2M /Tl} at half fiUing) in the ab- 
sence of impurities due to rare tunneling of the conden- 
sate phase introduced by the reservoirs, but can grow up 
to ^ M'^/Ep [Ep is Fermi energy in the wire) in the 
range of applicability of this model as the backscatter- 
ing increases. We impose spin dependent voltages in this 
model to examine the spin-charge separation by calculat- 
ing the average charge and spin currents and their zero 
frequency correlators and cross-correlator at temperature 
T . Each average current depends only on either charge 
Vc — or spin Vg = voltage combina- 

tion, respectively. At the voltages or temperature above 
the crossover both currents display a very different be- 
havior only weakly affected by the one electron scatter- 
ing. In particular, the spin voltage \Vs\ > F indeed re- 
verses direction of one of the spin polarized currents even 
for \Vc\ ^ \Vs\- The currents cross-correlator asymptoti- 
cally vanishes with increase of min |Vc,s|/F or T/F, while 
both diagonal correlators remain finite. This shows that 
the spin-charge separation in transport through the MHI 
occurs already in the Kondo regime at the energies less 
than Tl[< M). Below the crossover the spin polarized 
combinations of the average currents and their correlators 
approach in the linear voltage limit the correspondent de- 
coupled one-electron expressions [1| specified by the spin 
dependent voltages and a single parameter of the one- 
electron transmission. The latter is equal to 1/ (1-|-Fs/Fc) 
due to exponential enhancement of the one-electron re- 
flection T a/Tc by the charge condensate tunneling. 

Transport through a one-channel wire adiabatically at- 
tached to the left and right FL reservoirs is modeled by 
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a ID system of electrons, whose pairwise interaction is 
local and switched on by function (p{x) = 8{x)d{L — x) 
inside the wire of length L. Applying bosonization we 
can describe the charge and spin density fluctuations 
Pb{x, t) = ebdx4>b{x, t)/('\/27r), b — c,s, respectively, with 
(charge and spin) bosonic fields 0c, s- Without impuri- 
ties their Lagrangian symmetrical under the spin rotation 
reads 

- EpUb(pix)cos{2fj,bx/vb + V2(j)b{t,x))/{T:vF)] , (1) 

where vf{Ef) denotes the Fermi velocity (energy) in the 
channel. The parameter fic = fJ- varies the chemical po- 
tential inside the wire from its zero value at half-filling 
and = 0. The constants of the forward scattering 
differ inside the wire gb{x) = gb for x G [0, L] from 
those in the leads gb{x) = 1, and an Umklapp scat- 
tering (backscattering) of the strength Uc{Us) is intro- 
duced inside the wire. The velocities Vc,six) change 
from vp outside the wire to some constants Wc,s inside 
it. We can eliminate them rescaling the spacial coordi- 
nate Xoid in the charge and spin Lagrangians of ([T]) into 
Xnew = Jq dy/vc,s{y)- As a result, the new coordinate 
will have an inverse energy dimension and the length of 
the wire becomes different for the charge mode L — >■ 1 /T^ 
and spin mode L — > l/T'j^. Applying renormalization- 
group results of the uniform sin-Gordon model at en- 
ergies larger than Tl or we come to renormalized val- 
ues of the parameters in ([T]) . For the repulsive interaction 
when initially ffs > 1 > 5c: the constant Us of backscat- 
tering fiows to zero and g^ to 1, bringing the spin mode 
into the regime of the free TLL . The constant Uc of 
Umklapp process increases, reaching vp/vc at the energy 
cut-off corresponding to the mass of the soliton M if the 
chemical potential /i is less than M . Meanwhile, gc fiows 
to its free fermion value gc — 1/2. 

The spin-charge separation in Lagrangian ([T]) can be 
broken by additional one-electron scattering which en- 
tangles the spin and charge modes. We account for such 
a process by including a weak backscattering on a point 
impurity potential inside the wire < < L: 

2EpVimp (j)c{t,Xc) . (f>sit,Xs). , , 

— -cos( ^ +yo)cos( — ) > (2) 

where Xc.s ~ xq/v^^s, fo = 'fi + fJ-Xc incorporates a phase 
of the scatterer ip. The amplitude Vimp of the impurity 
potential determines transmission coefficient as 1/(1-1- 

1/2 ). 

imp/ 

Low energy model - An effective model for energies 
lower than some cut-off D' specified below has been de- 
rived 01 from the expression for the partition function 
Z associated to the combined Lagrangian ([T]) and ([5]) 
following Schmid Tc?]. Without impurities the spin and 



charge modes are decoupled. After integrating out (j)c in 
the reservoirs the charge mode contribution into Z de- 
scribes rare tunneling between neighbor degenerate vacua 
of the massive charge mode in the wire characterized by 
the quantized values of ^/2(f>c{T, x) + 2px = 27rm, m is 
integer. Variation of m by a — ±1 relates to passage 
of a (anti)soliton through the wire ((anti)-instanton in 
imaginary time r). The tunneling amplitude may be 
found as Pe"'*"/^^, sq = ^/A'P - M^, (a* < M) by map- 
ping [5||onto a free fermionic model or instanton tech- 
niques 171]. The latter also evaluates the pre-factor P = 
Cx Vi5'(si]rL/M2)i/4 with the constant C of the order of 
1. The parameter D' is a high-energy cut-off to the long- 
time asymptotics of the kink-kink interaction: F{t) = 
\iv{^Jt'^ + 1/_D'2} created by the reservoirs. It varies with 
/i and was estimated from the time scale of the instanton 
as D' ~ y/MTE at M = and D' ~ {M/fi)TL ff /i > Tp. 
A crucial modification to this consideration produced by 
the impurity under the assumption EpVimp ^ M en- 
sues from the shift of the m-vacuum. Since it is equal 
to {-l)"'EFAcos{<f)s{T,Xs)/V2), A = 2Vi„p cos .^/tt the 
neighbor vacua become non-degenerate. This can be 
accounted for by an auxiliary pseudospin variable with 
the correspondent Pauli matrix (T3. The energy split- 
ting becomes an operator a^^AEp cos{4>s{j,xo) / act- 
ing on the pseudospin, and every (anti-)instanton tun- 
neling reverses the as-value with the Pauli matrix cti. 
The interaction F{t) coincides with the pair correlator 
of some bosonic field ^c, whose evolution is ruled by 
the free TLL Lagrangian i2o[^c] ( uniform Lagrangian 
([!]) with no interaction). Then, by ascribing factors 
exp{^i6c{Tj,Q)/y/2) to the (anti-)instanton at the mo- 
ment Tj the functional integral for the partition func- 
tion is reduced Q to a standard Hamiltonian form Z = 
est X rr{e-^/'^} with 

n^nQ[<i3s{x)] +nQ[ec{x)]~ AEp€f3cos{(i)s{xs)/V2) 

- 2^6"^°/^^ 0-1 cos(6'c(0)/V2) . (3) 

Here (j)s{x) and 9c{x) are Schrodinger's bosonic oper- 
ators related to the variables (/)s(t, x) and 9c{t,x) of 
the functional integration. The free TLL Hamilto- 
nian Holcjisix)] {V-oidcix)]) is a function of the field 
(f)s{x) {9c{x)) and its conjugated corresponding to the 
free TLL action £o['/'s]('Co[6'c]), respectively. The model 
Q is equivalent to the initial one ([TJ at low energy. It re- 
lates to a Point Scatterer with internal degree of freedom 
in TLL and is solved exactly through fermionization. 

Fermionization - The Pauli matrices can be written as 
(Tq — (— 1)"+-^! J2p 7 ^"'^'''ipS.-y ■^ith Majorana fermions 
^1,2,3 and antisymmetrical tensor e : e^^^ — 1. Since the 
interaction in ([3]) is point-like localized and its evolution 
involves only the appropriate time-dependent correlators, 
we can fermionize it making use of: 

fry -»c(o) lEc' ■<t's(xs) 

MO) = -^J^^3e^,M0) = ^J^^ie—- (4) 
V Zn V ^TT 
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Here ipc,s{0) is the x = value of the charge (spin) 
fermionic field, respectively. These fields have linear 
dispersions taken after the related bosonic fields with 
momentum cut-offs (equal to the energy ones) D' and 
Ep, respectively. All states of negative energies are 
filled. Substitution of these fields into @ produces a 
free-electron Hamiltonian: 

a—c.s 

(5) 

where the interaction reduces to tunneling between the 
-tjjc.s fermions and the Majorana one ^ = ^2 ■ Here 
Ts = ^^{cos{ip)Vimp)'^ is the rate of the one-electron 
backscattering and the rate of the instanton tunneling is 

Application of voltages between the left and right 
reservoirs due to the shift of their chemical potentials, 
in general, different for electrons of opposite spin polar- 
izations cr = ± makes the system non-equilibrium and 
can be described with a gauge transformation (j)c.s ~^ 
<j)c,s — V2Vc,st in the real-time Lagrangian ([T|). Each in- 
stanton tunneling increases the condensate phase inside 
the wire and adds charge A(j)c/{\/2n) = 1 to the left 
reservoir. The correspondent change in the energy of the 
system equal to Vc causes a shift 6c/ V2 — >■ 9c/ V^ — Vet in 
the cos argument in Eq.Q. Both transformations then 
can be accounted for with the shifts of the charge and spin 
fermion chemical potentials in Eq. Q by Vc^s, respec- 
tively. Assuming below that both voltages are applied 
antisymmetrically and small enough, \Vc^s\ < Tj^ < M, 
we neglect their effect on the other parameters in the 
fermionized Hamiltonian ([3]). 

To find the charge and spin currents fiowing through 
the channel we notice that each antiinstanton realized 
by ^l^c{0) transfers charge 1 to the right reservoir with 
no transfer of spin. Then the charge current operator is 

Jc = -i[J dx^c^+,%F] = ^^/^c[^+{G) + i^cim = Jc- 

On the other hand, free passage of the spin current Js 
through the channel is affected by backscattering due to 
the spin field interaction in Hamiltonian "H in Eq. ([3|). 
This makes the spin current equal to Js = Vs/tt — js, 
where the backscattered spin current operator can be 
found as js = -V2SH/S^s = iVT's[i^+iO) + Vs(0)]e 

A crucial feature of the Hamiltonian "Hj? in Eq. ([5]) 
is that its interaction and the currents it creates con- 
tain two different Majorana components ?7a,=p(0) of each 
fermionic field ipa, respectively. These components are 
defined by V'a(0,t) — {r]a,+ — iT]a,-)/V2. Being orthog- 
onal at the same time both the 770, ±(i) components lose 
this property if taken at different times due to the applied 
voltages. Still, in Keldysh technique we need to use in 
non-equilibrium calculations both retarded and advanced 
cross-diagonal Green functions y^!^. (t) of free Majorana 
fermions vanish. Only the cross-diagonal Green func- 
tions g^'^_{t) = —g^'^^{t) are non-zero and equal to 



5>+1-H = [/((c.-K)/r)-/((a; + K)/r)]/2inthefre- 
quency representation, where / is the Fermi-distribution 
function. Then, from the Dyson equation the total cross- 
diagonal Green function 

G>1<^ = *V^(.9>+< -G^x + 9^+.a^G>;<) (6) 

reduces to the first product on the right-hand side, if the 
index X denotes ^ as the second field, or vanishes at all, 
ax = b+. The diagonal total Green function G^_^_ 
coincides with the free one. 

Average currents - As follows from Eq. ^ the aver- 
age of the operator < ja > can be written as < ja >— 
—i2Ta J dojg^ ^ _{oj)G^{uj)/ {2Tr) and depends only on 
the corespondent voltage Va, since the advanced Green 
function does not contain information about voltages. 

It is related to the free Green function 5^^^'' — 2/ (w^iO) 
through the correspondent Dyson equation with the self- 
energy S^^^-* = zLiT. Substitution of the expressions for 
both Green functions results in 

^ . ^ 2r J f /(^)-/(^) 

<Ja>= duj 2 I /ir2 ' 

The average charge and spin currents Jc,s are equal to 

< jc > and Vs/n— < js >, respectively. At low tem- 
perature T < r the integral in Eq. ([7]) converges to 

< (2ra/7r)arctan(K_/(2r))._ Theii the aver- 
age spin polarized currents = {Jc + crJs)/2 below 
the crossover F approach the one-electron expressions 
Ja = DVa/{2TT), where the spin independent transmit- 
tance D = ^ demonstrates renormalization of the ini- 
tial amplitude Vimp in Eq. ^ into ^F^/Fc by the 
interaction inside the MHL Above the crossover the 
charge current saturates at Fc while the spin current 
grows up as Vs/n — Tg- Therefore, for \Vc\ » \Vs\ 
one of the conductances Ja/Va becomes negative as 
\^s\ ^ ''■r suggesting emergence of the spin-charge sep- 
aration. Similarly, the separation may be expected at 
T :g> F. Indeed, the linear bias charge conductance is 
Gc = FcV''(l/2-hF/(7rT))/(7r2T), where ^p'{x) is the 
derivative of the di-gamma function, ip'{l/2) = tt'^/2. 
The high temperature asymptotics of both conductances 
Gc = Tc/(2T) and Gs ^ l/n - F J(2T) are defined by 
different F-parameters and independent of each other, 
and the condition 1 ^ = Gc/Gs on the current 
reversing voltages is satisfied. We further examine corre- 
lators between the charge and spin currents. 

Currents correlators - The zero-frequency current cor- 
relators 6"^ Jab are related to the current operators corre- 
lators 5'^ jab — J dtexp[—iujt] <ja{t)jb{0)>,uj — >• in the 
following way S^Jab — ^S'^jab, where ± stands for diago- 
nal and cross-correltors, respectively. The diagonal s pin 
current correlator at finite temperature also includes 

additional terms (5^ Jss =2T/7r-4r9v/ <js > +S'^jss- Ap- 
pearance of these terms may be easily understood recall- 
ing that the fluctuation-dissipation theorem claims non- 
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zero current fluctuations even in the absence of backscat- 
tering at non-zero temperature. 

Substituting expressions for the current operators ja^b 
into their correlator and then spUtting the correlator into 
pair-wise correlators of the Majorana fields by applying 
Wick's theorem we find that 



ab 



(2) 



(8) 



Then the cross-correlator of the two currents equal to 
-lil^ in Eq. ^ follows from Eq. © as 



6' Jr.. = 



4rj, 

27r 



rfc.[Gf(a.)]^5>+-(^)3s<+-(^) 



n..e. [/(^)-/(^)] 



{oo-2iTf ■ 

It is vanishing with increase of both voltages or temper- 
ature under the integral in Eq. ^ due to the analytical 
structure of the advanced Green function. In particular, 
in the limit of |V"c,s|/T ^ 1 it takes the form 

4r,r, sign(V^,T4)F 



TT [4r2 + 



V — min 

a—c,s 



Va\ 



(10) 



which approaches zero oc TcTg/V as V becomes much 
larger than 2r. In general, the integral in Eq. (jH]) is ex- 
pressed in terms of derivatives of the di-gamma function. 
This expression shows that at high temperature T ^ 2r 
the cross-correlator is vanishing as S^Jcs = li^^^H^/'^) + 
^P^^\l/2)\TcT,VcVs/il2Tr^T^) if 1 » maxl^e.^l/T. 

The diagonal correlators of the charge and spin cur- 
rents also include the second term in Eq. ([8]). To 
find it we notice that the Green function of the 
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single Majorana operator reduces to G^ 
with the self energy I]> = -iEa.±^af{{~uj ± Va)/T) 
and g>++(-cj) = -«E± /((^ ± Substitution of 

these expressions into Eq. ^ gives us the diagonal zero- 
frequency correlator of the charge current in the following 
form 



2r, 



T ' 



/(- 



T 



/(- 



-UJ - 



4r2 



T 
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and the similar expression for correlator of the spin 
backscattered current 6^jss by exchanging the "c" and 
"s" -indexes. At zero temperature the latter correlator 
coincides with 6^Jss and hence both shot noises of the 
charge and spin current (a = c, s' 
follows: 



rise from Eq. ([Tl]) as 



2rj, 



X arctan 



/ max{|T4|}j,^c 

^ 2r 



4r^lK| 



4K' + 4r2) 



(12) 



In the low voltage limit both expressions coincide S^Jcc = 
S'^J,^ = {l/Tr)D{l - D)TJia^{\Vb\}b=c,s- Combined with 
the low voltage cross-correlator in Eq. (|10p they give 
shot noise of the two spin polarized currents [8] in the 
one-electron form S^J^a' = {Saa' /t^)D{^ — ^)|Kt|- 

As both voltages increase above the crossover the shot 
noise in the two currents grows and saturates at two dif- 
ferent values S^Jaa = Ta- Together with demonstrated 
before suppression of the cross-correlator of the cur- 
rents this confirms the spin-charge separation above the 
crossover but yet below Tj^ . Similarly, the charge current 
noise d^J^c = + 2iP^^\l/2)T/{Tr^T)] remains finite 
at high-temperature, while the spin noise grows linearly 
S^Jss = 2T/7r-Ts = 2TGs. 

Toulouse limit in Kondo dot transport - This model has 
been written Q as a formal generalization of the physi- 
cal model of Kondo dot transport. It describes tunneling 
between two branches, a = r,l, of ID chiral fermions 
ipa,cT{x) carrying spin a and propagating in the right and 
left reservoirs, respectively, with the tunneling Hamilto- 
nian 



.^h—r.l (t,(t' ,7 



/;-''f^<,(ox^,V6,.'(o) 



(13) 



where are the Pauli matrices for the Kondo dot spin. 
The choice of the parameters Jf" = 2ti, J^^^ = 0, J^;,^ = 

J"'*" corresponds to the Toulouse limit model solvable 
through bosonization and Emery-Kivelson rotation. Un- 
der additional restriction J^'' = — J^' its Hamiltonian 
takes the form of H in Eq. ([3|) with the field (j>s sub- 
stituted by its dual 6s, which describes tunneling of the 
unit of spin instead of its backscattering, and further co- 
incides with Hp in Eq. ([5]) after the fermionizaton with 
r, - EF[r/]y i87r),rs = i?F[J7]V(87r). Both oper- 
ators jc,s now describe tunneling of the charge and the 
spin. Interchanging the direct and backscattered spin 
currents in the above calculation we apply their results 
to the Toulouse limit model. In particular, the low en- 
ergy linear voltage dependence of the average charge and 
spin currents is defined by different transmittances equal 
to D and 1 — D, respectively. This feature can not be 
accounted for in the one-electron transport model. In- 
deed, in this model the spin cross-diagonal transmission 
and reflection coefficients become zt{D — l /2), and one of 
them would be negative unless D ^ 1/2. This indicates 
independent tunneling of spinous and break of the low 
energy electronic PL description. 

In summary, we have shown that in spite of influence 
of the PL reservoirs and a weak one-electron scattering 
the essential properties of the spin-charge separation in 
transport through a MHI of finite length appear if either 
T or both are above a crossover energy (^ Tl), be- 
low which, however, the transport becomes one-electron. 
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